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Abstract 

A class of predictive SO (10) grand unified theories with highly asymmetric mass matrices, 
known as lopsided textures, has been developed to accommodate the observed mixing in the 
neutrino sector. The model class effectively determines the rate for charged lepton flavour vi- 
olation, and in particular the branching ratio for fi 67, assuming that the supersymmetric 
GUT breaks directly to the constrained minimal supersymmetric standard model (CMSSM). 
We find that in light of the combined constraints on the CMSSM parameters from direct 
searches and from the WMAP satellite observations, the resulting predicted rate for /i — 67 
in this model class can be within the current experimental bounds for low tan/3, but that 
the next generation of yU — > 67 experiments would effectively rule out this model class if LFV 
is not detected. 
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1 Introduction 



Neutrinos have been observed to oscillate between flavour states P]-|HI; which implies neu- 
trino mass and mixing. In addition, the combined observations suggest that both the at- 
mospheric and solar mixing angles are nearly maximal, known as the large angle mixing 
solution (LMA). Interestingly, the LMA solution implies that the lepton mixing scenario is 
radically different from the quark sector. From the low energy point of view, we should 
expect that the neutrino mass inducing dimension 5 operator (HHLL) would be the first 
observable signal beyond the renormalizable operators that compose the standard model. 
Furthermore, the smallness of the inferred masses suggests that the mechanism responsible 
for neutrino mass is distinct from purely electroweak physics, and could naturally arise from 
physics at a very high scale. The dimension 5 operator can be induced by adding three 
heavy gauge singlet Majorana fermions (one for each generation) to the standard model. 
Upon integrating the heavy Majorana fermions out at their associated scale, small neutrino 
masses are induced after electroweak symmetry breaking. This neutrino mass generating 
technique is the see-saw mechanism jUj. 

While the see-saw mechanism is an economical and natural way to understand the small- 
ness of the inferred neutrino masses, there are many possible methods of implementing it, 
and therefore detailed neutrino observations can be used to constrain GUT models. Perhaps 
the most elegant GUT uses the grand unifying group SO(IO) in four spacetime dimensions. 
The spinor representation of SO(IO) is 16 dimensional, which accommodates all the helicity 
states of one fermion family plus an extra singlet degree of freedom for a Majorana neutrino. 
The generations are simply three copies of the spinor representation. Since GUTs relate 
quark and lepton masses and mixings, it is perplexing from a model building perspective as 
to why lepton mixing is so different from that in the quark sector. More specifically, it is of 
interest to understand why jf/^sl of the MNS matrix is so much larger than \ Vcb\ of the CKM 
matrix. Over the last few years a number of models have been developed to address this 
difference [IHl^lini- Recently, a particularly interesting and highly successful class of super- 
symmetric SO (10) GUTs has emerged that makes use of asymmetric mass matrices known 
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as lopsided textures jTUl HH [121 • these models, the charged lepton sector is responsible for 
the large atmospheric mixing angle while the Majorana singlet neutrino matrix has a simple 
form that results in the large solar mixing angle. Throughout this paper we will refer to 
these models as the AB model class |lUj . 

After GUT breaking, these models reduce to the R-parity conserving minimal supersym- 
metric standard model (MSSM) with specific model dependent relationships amongst the 
Yukawa couplings. In addition to the constraints already provided by the neutrino physics 
(and the demand that these models reproduce all the low energy physics of the standard 
model), the WMAP satellite observations ^U] provide strong constraints on the available 
supersymmetric parameter space if the lightest supersymmetic particle (LSP) is assumed 
to compose the dark matter fl^ 121] ■ Assuming the constraints on the CMSSM from 
the WMAP data, the definite flavour structure of the AB models will result in specific soft 
supersymmetry breaking parameters. Therefore, the AB model class gives well defined pre- 
dictions for lepton flavour violation and in particular fi 67. It is of considerable interest 
to determine how the lepton flavour changing neutral current bounds restrict the CMSSM 
parameters for the AB model class in light of the WMAP data. 

We organize this paper as follows. In section |21 we outline the essential details of the AB 
models, the supersymmetric parameter space, and the calculation for fi 67. We consider 
fi —>■ ej since at the present time, with the current bound [25^ of BR(/i — > 67) < 1.2 x 10~^^, 
this process gives the strongest constraints on lepton flavour violation in the class of models 
that we discuss. Furthermore, the MEG experiment at PSI [2^] expects to improve on this 
bound with the expected sensitivity of BR(/x — >■ 67) ^ 5 x 10^^^. This experiment will 
provide stringent limits on models with charged lepton flavour violation. In section |21 we 
display our numerical results with the combined constraints from /i — > 67 and the WMAP 
satellite observations, and in section |3] we present our conclusions. The appendix provides 
further calculational details. 
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2 The AB Model Definition 




e 



The AB model class is based on an SO(IO) GUT with a U(l)xZ2xZ2 flavour symmetry and 
uses a minimum set of Higgs flelds to solve the doublet-triplet splitting problem [TUl HH [121 • 
The interesting feature of these models is the use of a lopsided texture. The approximate 
form of the charged lepton and the down quark mass matrix in these models is given by 

/nnn\ /oOO^ 

Yd ~ (T . (1) 
a 1 J yOely 

where a ~ 1 and e ^ 1. As pointed out by the authors of [TU], this asymmetric structure 
naturally occurs within a minimal SU(5) GUT where the Yukawa interaction for the down 
quarks and leptons is of the form AjjSjlOjSH (5h denotes the Higgs scalars). In an SU(5) 
GUT, the left-handed leptons and the charge conjugate right-handed down quarks belong to 
the 5 while the 10 contains the charge conjugate right-handed leptons and the left-handed 
down quarks. Therefore the lepton and down quark mass matrices are related to each other 
by a left-right transpose. Since SU(5) is a subgroup of SO (10), this feature is retained in an 
SO(IO) GUT. This lopsided texture has the ability to explain why \U^s\ >> \Vcb\- Making 
use of this observation, the AB models contain the Dirac matrices U, N, D, L for the up-like 
quarks, Dirac neutrino interaction, down-like quarks, and the leptons respectively 



u 



V 














e/3 





-e/3 


1 



Mu, N 



^ ri ^ 

-e 

e 1 y 



(2) 



D 






6 


6 







-e/3 



a + e/3 
1 



D, 






6 


6 


-e 




cr + e 1 



(3) 



3 



where 

Mu ^ 113 GeV, Md ~ 1 GeV, 

a = 1.78, e = 0.145, 

(4) 

6 = 8.6 X 10-3, S' = 7.9 X lO^^, 

= 126°, T] = 8x 10'^. 

Dimensionless Yukawa couplings Yu,Yn,Yd, and Ye can be extracted from the Dirac 
matrices. The given values of Mo and Mjj best fit the low energy data with tan/3 ~ 5. It 
should be noted that larger values of tan (3 are easily accommodated by altering the values 
of M[7 and while retaining accurate fits to the low energy data after renormalization 
group running. The lopsided texture of the AB model class nicely fits the large atmospheric 
mixing angle; however, in order to obtain the large solar mixing angle a specific hierarchical 
form of the heavy Majorana singlet neutrino matrix needs to be chosen [TH I12j. namely. 



^ h'^rf —ber] ar] ^ 



An. (5) 



—bet] — f 

y ar] — e 1 J 

where the parameters e and t] are as defined in equation The parameters a and b are of 
order 1 and An ~ 2 x 10^^ GeV. Since the Majorana singlet neutrino matrix is not related 
to the Dirac Yukawa structure, it is not surprising that this matrix should take on a form 
independent from the rest of the model. Once these choices have been made, the AB model 
class is highly predictive and accurately fits all the low energy standard model physics and 
the neutrino mixing observations. 

It should be emphasized that all these relations are defined at the GUT scale and are 
therefore subject to renormalization group running jOj. If we conservatively assume that the 
GUT symmetry breaks directly to the standard model gauge symmetries, SU(3) xSU(2)xU(l), 
and that supersymmetry is broken super-gravitionally through a hidden sector in a flavour 
independent manner, the AB model class will give well defined predictions for charged lepton 
flavour violation. There may also be significant contributions to the off-diagonal elements 
from renormalization group running between the GUT and gravity scales [2HI ■ Since the 
particulars of GUT and supersymmetry breaking - as well as the possibility of new physics 



above the GUT scale - can have model dependent effects on the branching ratio for /i 67, 
we do not consider an interval of running between the GUT and gravity scales. 

The specific model predictions for the Dirac Yukawa couplings and the form of the Ma- 
jorana singlet neutrino matrix will feed into the soft supersymmetry breaking slepton mass 
terms through renormalization group running, generating off diagonal elements that will con- 
tribute to flavour changing neutral currents . The amount of flavour violation contained 
in the AB model class can be examined through the branching ratio of the process n 67. 



3 Numerical Results for ^ 

After GUT and supersymmetry breaking, we have the constrained minimal supersymmetric 
standard model (CMSSM) with heavy gauge singlet neutrinos to make use of the see-saw 
mechanism. The leptonic part of the superpotential is 

W = e„^i/d"EYEL^ + e„/jKNYNL^ + ^NMnN (6) 

where Ye, Yn are Yukawa matrices, and Mn is the singlet Majorana neutrino mass matrix. 
The totally antisymmetric symbol is deflned ei2 = +1. We explain our notation in detail in 
the appendix. On integrating out the heavy singlet neutrinos, equation (0) reduces to 

W = e,^i/d"EYEL^ - ^z/^m^i. (7) 

where 

= yY^MN'YNsin2/3 (8) 

is the see-saw induced light neutrino mass matrix. The coefficients (3 and v are deflned in 
terms of Higgs flelds expectation values by 

^ = {Hlf + {H'^y = (174 GeV)^ , tan/3 = |||. (9) 

The neutrino mass matrix, equation (jH)), is in general not diagonal and this is the source of 
lepton flavour violating interactions. 
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We assume that supersymmetry is broken softly in that breaking occurs through operators 
of mass dimension 2 and 3. The soft supersymmetry breaking Lagrangian relevant to LFV 
studies is 

(^-Be^^H^Hi - ^NB^N + c. c.^ 

(^-e^aH^EAEL^ - eapH^NA^t'' + c. c.) 
--MiBB - ^AW^VF" + c. c. ) (10) 

(see the appendix for the notational details). The CMSSM assumes universal soft super- 
symmetry breaking parameters at the supersymmetry breaking scale, which we take to be 
of order the GUT scale, leading to the following GUT relations: 

= "^1 = "^1 = "^0 • I> (11) 

"^Hd = "^Hu = "^o> (12) 

Ae = An = 0, (13) 

Ml = M2 = my 2 (14) 

where mo and mi/2 denote the universal scalar mass and the universal gaugino mass respec- 
tively (I is the 3x3 unit matrix). We conservatively assume that the trilinear terms Ae and 
An vanish at the supersymmetry breaking scale. 

We run the parameters of the CMSSM using the renormalization group equations (see 
appendix) working in a basis where the Majorana neutrino singlet matrix is diagonal, in- 
tegrating out each heavy neutrino singlet at its associated scale. After integrating down 
to the electroweak scale, we rotate the Yukawa couplings to the mass eigenbasis. In order 
to understand the origin of flavour violation in this model class, wc first give a qualitative 
estimate. The leading log approximation of the off-diagonal slepton mass term is given by 

(A-a., « - J^-^SIY.'YJ In (^^) , (15) 



(assuming that the trihnears vanish at the GUT scale), and using this approximation together 
with mass insertion techniques [221 ISO] , the branching ratio for /z — > 67 is 

2 
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87r2 ^ An 

where rris is a typical sparticle mass. We see that since the flavour structure of the AB 
model class is specifled so precisely, the branching ratio for — 67 is well determined. 
In our calculation of the decay rate, we use the full one-loop expressions derived from the 
diagrams in flgure[T](see the appendix for more details). 

The WMAP satellite observations 21^ combined with constraints from b ^ s'y and LEP 
direct searches |31j strongly limit the available CMSSM parameter space if the LSP composes 
the dark matter [^1^ 1^. In addition to these constraints, realistic supersymmetric GUT 
models must also survive LFV bounds, such as the limit on /i — ej. In particular, using all 
of the available bounds, both cosmological and laboratory, we can further restrict the AB 
model class. 
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Figure 1: Feynman diagrams contributing to /i — 67. 
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In figure 121 we show contours of tlie brancliing ratio /i — 67 in tlie mi/2-mQ plane for a 
variety of tan (3 with the /i parameter both positive and negative. The parameters of the AB 
model class have been chosen such that all the low energy predictions fit the standard model 
data, and we have chosen a = 1 and 6 = 2 for the Majorana singlet neutrino mass matrix 
given in equation (0). As indicated in there are a number of possible model choices 
for the Majorana singlet parameters a and b that are consistent with the LMA solution. 
However, we find that the rate for — 67 is largely unaffected by the allowed range J2] 
for these parameters. Panel (a) demonstrates the lepton flavour bounds for tan (3 = 5 with 
/i > . The small line-like shaded area in the lower part of the panel is the allowed region 
from the combined WMAP and laboratory limits. The remaining panels show that the 
contours of constant branching ratio migrate to the right of the plots (i.e. to high values of 
mi/2 and mg) as tan/5 is increased. In each case we overlay the approximate WMAP and 
laboratory constraint bounds represented by a shaded region P2j. The choice for the sign of 
/i is indicated in each panel. As tan (3 is pushed up, larger portions of the parameter space 
become excluded. This is an expected feature since the branching ratio is proportional to 
tan^ (3. Notice that by tan /3 ~ 25, /i > 0, the branching ratio allowed contours no longer have 
a significant overlap with the WMAP region. As a result, we find that the AB model class 
is consistent with the current experimental bound on /i — > 67 for low tan/5 (i.e. tan/3 < 20) 
for /i > 0. For completeness, in panels (b) and (e), we show two cases where /i < 0. The 
branching ratio of /i — 67 is largely insensitive to the sign of fi, however the WMAP region 
is moderately affected |21] . A small part of the allowed WMAP region is currently permitted 
for larger tan/3 (i.e. ~ 35) as indicated in panel (e). The upcoming limits ^jj that MEG 
will establish, BR(/i 67) ^ 5 x 10^^^, will effectively rule out this model class if LFV is 
not seen. Interestingly, if LFV is seen at MEG, this model will suggest that tan/3 is low 
based on flavour bounds alone. 
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Figure 2: Contour Plots of BR(/i — >• ej) in the mo — mi/2 plane: Panels (a),(c),(d), and (f) 
show the contours of the branching ratio for tan f3 = 5, 15, 25, 50 respectively with fi > 0. 
Panels (b) and (e) show the contours with tan f3 = 10, 35 respectively with fi < 0. In all 
cases the shaded region corresponds to the approximate combined WMAP and laboratory 
constraints. 

4 Conclusions 

The AB model class [TUl ITU IT^ . based on a U(l)xZ2xZ2 flavour symmetry, is a highly 
successful and predictive GUT scenario. This model class has the ability to accommodate 
all the observed neutrino phenomena and reproduce the low energy physics of the standard 
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model. If it is assumed that supersymmetry is broken via mSUGRA and that the GUT 
breaks directly to the CMSSM, the AB model class is highly restrictive and hence allows for 
a precise determination for the rate of charged lepton flavour violation. In particular, we 
examined the process fi — > ej, since at the present time this flavour violating muon decay 
channel gives the strongest constraints on flavour changing neutral currents in the lepton 
sector. 

As the WMAP satellite data [21] and laboratory direct searches have already severely 
restricted the available CMSSM parameter space, the /x — > 67 flavour bounds allow a strong 
test of the AB model class. We flnd that given the current bounds on /i — 67, BR(yU — * 
67) < 1.2 X 10^^^, the AB model class favours low tan/5 (i.e. < 20) with /i > 0, however, 
there is a small region that is not excluded for tan/3 < 35 with the sign of fi negative. If 
MEG at PSI j2n] does not detect a positive LFV signal, BR(/i ^ 67) < 5 x 10-^^ the 
AB model class will be effectively ruled out, given our conservative assumptions concerning 
GUT and supersymmetry breaking. It remains an open question as to whether or not other 
supersymmetry and/or GUT breaking schemes within the AB model class will be able to 
avoid these flavour violating bounds. 
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6 Appendix 

In this section we wish to clarify some of the calculational details. We carefully establish our 
notation and conventions. Also, we include the full one loop amplitude for the rate fi e'j 
that we used in our calculations. Formulas similar to those given in subsections I6.2|IU?^ can 
be found in 
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We express the supersymmetric Lagrangian using the 2-component Weyl formahsm. — 
(L", L2 , L'^)'^ denotes a column vector in generation space containing the SU(2) doublet 
lepton chiral supcrfields; 1,2,3 are generation labels, and a = 1, 2 arc the SU(2) indices. E = 
(£^1, E2, Es) denotes a row vector in generation space containing SU(2) singlet charged lepton 
superfields. The gauge singlet neutrino chiral superfields are denoted by N = (A^i, N2, N^). 
Similarly, for the quark superfields: = (Q^, Q^, Q^)'^ denotes the SU(2) doublet, = 
u = {ui,U2,U3)^ , = d = (di, ci2, (is)^; and the SU(2) singlet quark superfields are 
U = {Ui,U2,U3), D = {Di,D2,D3). H^, are the SU(2) Higgs doublet superfields of 
opposite hypercharge with the standard components: H^^^ — H^, H^^^ = H^, H^^^ — H^, 
= H^. The corresponding scalar components of the superfields are written respectively 
as L", — Jj'^ — e; E; N; Q", Q-*^ = u, = d; D; U; (all are vectors in generation 
space). The fermionic components of the Higgs superfield, the Higgsinos, are denoted as H^, 
H^. The superpotential W is given by 

+ iie^pH^H^ (17) 

where Yg, Yn, Yd, Yy are Yukawa matrices Mn is the singlet Majorana neutrino mass 
matrix, is the Higgs parameter that breaks the U(l) Pecci-Quinn symmetry, and the totally 
antisymmetric symbol is defined 612 = +1. The soft supersymmetry breaking Lagrangian is 

>Cbreaicing = -(5«^L"W?L^ - Em|Et - Nm|Nt 

(^-be^f^H^H^ - ^NBj^N + c. c.^ 

(-e«^i7d"EAEL^ - e^^pH^NA^l^ + c. c.) 
(^-eapH^bAuQf' - SapH^lJA^Q'' + c. c.) 

(--MiBB - -Msiy^iy" - -MaG^G^ + c. c. ) (18) 

11 



where B denotes electroweak U(l) gaugino field; W"", a = 1,2, 3, denote electroweak SU(2) 
gaugino fields; &, b = 1, ...,8, denote strong interaction, SU(3), gaugino fields; m?, m|, 
m|, m|, m|, m|, B^,, Ae, An, Ad, Au, m^^, m^^, b. Mi, M2, M3 are the supersymmetry 
breaking parameters, and at the GUT scale: 

2222222t /inN 
"^L = "^E = "^N = "^Q = "^D = "^U = "^0 ■ I. (19) 

"^L = "^L = (20) 

Ae = An = Ad = Au = 0, (21) 
Ml = M2 = M3 = mi/2 (22) 

where mo and mi/2 denote the universal scalar mass and the universal gaugino mass respec- 
tively (I is the 3x3 unit matrix). After running the CMSSM RGEs (see subsection 16. 6p . we 
rotate all the Yukawa couplings to a diagonal basis, and in particular the lepton sector. 

Ye ^ U^YeVJ, = diagonal, (23) 

m? VEm?vL (24) 

m| ^ U*m|U^, (25) 

Ae ^ U*AeV^. (26) 

Not all of the bi-unitary rotation matrices can be absorbed away through the field re- 
definitions as the left-handed neutrinos become massive below the see-saw scale and after 
electroweak symmetry breaking. 

6.1 fi parameter 

The scalar potential of the Higgs fields is given at its minimum by 

+ b{H',){H'^)+b*{H',){H^,) 

+ '-^^{wy-wYy (27) 
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where gi, g2 are respectively U(l) and SU(2) gauge coupling constants. We can use the 
SU(2) gauge transformation freedom to choose the vacuum expectation value of the charged 
Higgs field {H^'^ = 0; then it follows that also {H^) = at the minimum of the Higgs 
potential. Therefore, we are left with only the neutral Higgs fields of equation (jTTj) . The 
conditions that the minimum of the potential V breaks the electroweak symmetry properly 
are 



/i^ + m^^j + b tan P 
+ m^^ + h cot (3 



--ml cos 2/3, 
-mg cos 2/? 



(28) 
(29) 



where mz is the mass of the Z-boson. After eliminating the terms containing h we obtain 
the tree level parameter relation, 

,2 ^2 i„„2 



2 1 2 "^Hd - "^H, tan^ (5 
^ = l-tan^/3 • 



(30) 



6.2 Neutralinos 

The neutralinos Xi? X2' Xs? X4 are mass eigenstates of the neutral gauginos B, and 
neutral Higgsinos if^, H^. The neutralino mass Lagrangian is given by 

/ B \ 
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B ] Mr 
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—niz cos P sin 
y mz sin /? sin 9w 
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mz cos /5 cos 6'w 
-mz sin /5 cos 6'w 



— mz cos (3 sin 6'w 



(31) 



mz sin (3 sin 9 
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mz cos (3 cos 6'w — w^z sin (3 cos 6'w 








(32) 
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An orthonormal rotation leads to the mass eigenstates: 

/ B \ 



Or 



(33) 



where One is a real, orthogonal matrix. The mass matrix (j32|) can therefore be decomposed 
in terms of real mass eigenvalues, M^o, a = 1, 2, 3, 4, 

= Oj.diag (m^o M^o M^o M^o) One, (34) 

and can be rewritten as 

^ = -\ilM&lxl (35) 



a=l 



6.3 Charginos 

The charginos are mass eigenstates of the charged SU(2) gauginos and charged Higgsinos, 

w- 



c. c. 



where 



w~- 



V2 ■ 



Mr 



and the mass matrix is 

M2 \/2mw cos (3 
v^mw sin P /i 
(mw is the VT-boson mass). The mass eigenstates are given by 



(36) 



(37) 



(38) 



Xi 
X2 



Or 



O 



R 



(39) 



where Or and Ol are real orthogonal matrices, and they can be chosen so that the mass 



eigenvalues M--, M-- are positive, and 



Mc = 0;^diag (M-- M--) Ol, 



(40) 
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Equation ()36|) can be written as 



(41) 



6.4 Sleptons 

Masses of the charged sleptons are given by the Lagrangian 



'RR-" 



(42) 



with the mass matrices 



m 



LL 



m 



RR 



m 



RL 



nil + "^z COS 2/5 l^sin — - 

+ m| — m| cos 2/5 sin^ 6'w ■ I, 

vcosP ^ 
— /iirii tan p H ^^Ae 



where 



V2 



mi = diag (m^ mi^ 



(43) 
(44) 
(45) 

(46) 



and m^, rrii^, are electron, muon, and tau masses respectively. The above Lagrangian 
written in terms of mass eigenstates /i, /e (six complex scalar fields) is 

6 



6=1 



with 



h 
h 
h 
h 

and Ujf is a complex unitary matrix defined by 



/ft ^ 

62 

v4v 



(4J 



2 2t 



m 



RL 



m 



RR 



U J diag [m]^ m| m| m| m| m|) Uf. 



(49) 
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Similarly, the light sneutrinos (the heavy singlet sneutrinos are ignored since they have 
decoupled well above the weak scale) 



(50) 



where 



M? = m? + -m| cos 2/3 • I. 



(51) 



The sneutrino mass Lagrangian written in terms of mass eigenstates ni, fi2, (three complex 
scalar fields) reads 

" ~ (52) 



£ = - ^ ml^hlhb 



6=1 



with the mass eigenstates defined by 



n2 

and Ufi is a complex unitary matrix satisfying 



M2 = utdiag {ml ml ml)V^. 



(53) 



(54) 



6.5 Lepton Flavour Violating Interactions 

The interactions leading to the lepton flavour violating process Ij —>■ h + ■y involve two 
effective Lagrangians: neutralino-lepton-slepton and chargino-lepton-sneutrino. Written in 
the mass eigenbasis they are 

^ = E E E ^l^admi + N^ineixi + c c. (55) 



and 



1=1 a=l 6=1 



3 2 3 



^ = E E E CLhE^Xa + &le^x'a + C C. 



(56) 



i=l a=l 6=1 
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where 



^1 = ^ (tan^w (U^);, (O.,),, + (U^);, (0.e)„2 - ^^^^^ (Uf);(,^3) (One).3) , (58) 
and 

CL = ^ (Ql).. (U,):, , (59) 

y/zniT^ cos p 

Cfa, = (Or), JU,),V (60) 

The on-shell amphtude for Ij — ^ /j + 7 can be written in the general form 

M = ee^h {p - q) [im^M^^q, (AlL + ArR)) (p) ; (61) 

here we have used Dirac spinors /j {p — g) and Ij (p) for the charged leptons i and j with 
momenta p — q and respectively; L = (1 — 7^) /2 and R = (1 + 7^) /2. Each of the 
dipole coefficients and Ar have contributions from the neutralino-lepton-slepton and the 
chargino-lepton-sneutrino interaction, namely, 

A^ = A^^^ + At\ (62) 

= + (63) 
where , A^^ , ^ , A^^^ can be evaluated from the Feynman diagrams in figure d 







327r2 ^ 
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1 / /ML\ M- /ML , , 

' ^L* T Xa \ I ^R* _Xa J Xa \ \ f f!Cr\ 
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(66) 
(67) 



The functions Ji (x), J2 (x), J3 (x), J4 (x) are defined as 



^ , , 1 — 6x + 3x^ + 2x^ — 6x^ Inx , , 

-^1(3;) = — -I , (68) 

6 (1 — x) 
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1- + 2x\nx 

J2[x) = -3 , (69) 

(1 -x) 



2 + 3x — 6a;^ + x'^ + 6x In a; 

6(1 -x)' 
-3 + 4x - x^ + 21nx 



■^3(3;) = ;:7^ ^4 ' (70) 



J4(x) = -3 . (71) 



Finally, the decay rate for L — > / ■ + 7 is given by 



r{l-^l-+l) = j^ml{\A,f + \A^f), (72) 
and i = 1, J = 2 for /i — i> e + 7. 

6.6 Renormalization group equations (RGEs) 

The general form of the supersymmetric renormalization group equations [221 E21 120] are 

dt 167r2 ^ ^ 

where X is any of gi, g2, gs, Yn, Ye, Yu, Yd, Mi, M2, M3, m^^, m|^, m?, m|, m|, m|, 
m?j, m|. An, Ae, Ay, Ad, and the dotted quantities are listed below: 

gi = ngl (74) 
g2 = gl (75) 
g3 = -3gl (76) 

(77) 

Yn = Yn (^-g^I - Sg^I + 3Tr (y{jYu) I + Tr (yJjYn) I + 3YJ,Yn + yJ^Ye) , (78) 
Ye = Ye (^-3gfl - 3gll + 3Tr (yJ^Yd) I + Tr (y^Ye) I + 3YJ,Ye + yJ,Yn) , (79) 

Yu = Yu(-y^??I-3(72l-^^2l + 3Tr(Y{jYu)l + Tr(Yj,YN)l 

+ 3Y{jYu + Yf,YD), (80) 
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Yd = YD(^--^?I-3(72l-y^3'l + 3Tr(Yj,YD)l + Tr(Yt,YE)l 
+ 3YiYD + YiYu), 

Ml = 22glM,, 

M2 = 2glM2, 
Ms = -GglMs, 

S^mi^- m^^ + Tr (m| - 2m| + m| - m? + m|) , 

m^^ = 6IV (m^YfjYu + Yfjm|Yu + m^^ YfjYu + A}; Au) 
+2Tr (m? yJjYn + Y^mlYN + Y^Yn + AJ,An) 
-2gfMi-6giMl + gfS, 

ml^ = 2Iv(m?Yj,YE + Y;,m|YE + m^XYE + At,AE) 

+6Tr (m^Yj^YD + Yj,m|YD + yJ^Yd + AJ^Ad) 
-2glM', - QglMl - glS, 

lii? = m?Y^YE + Y^YEm? +m?Yj,YN + YjjYNm? 
+2Yj,m|YE + 2m^^Yj,YE + 2AJ, Ae 
+2Yjjm|YN + 2m^^Yj,YN + 2A3,An 
-2glM!l-6glMil-glSI, 

m| = 2m|YNYjj + 2YNYjjm| + 4YNm? yJ, + 4m^^YNYjj + 4AnAJj, 
m| = 2m|YEYj, + 2YEY^m| + 4YEm?Y^ + 4m^^YEY^ + 4AEA^ 

m| = m|YiYu + Y;jYum| + 2Yim|Yu + 2m2^Y;jYu + 2A{jAu 

+ Y{)YDm| + 2Yfjm|YD + 2Tr4^ Y{,Yd + 2A{)Ad 
-^^i^M,^I - 6p,^M|l - fglMil + ^glSl, 
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ihl = 2m?jYuY{j + 2YuY{jm| + 4Yum|Yi + 4m^^YuY{j + 4AuAjj 

-fgMl - fgMl - lalSI, (92) 

Ihl = 2m|YDYj, + 2YDY]3m? +4YDm|Yj, + 4m^^YDY]3 + 4ADA|, 

-IgMl-fsMl + lglSl, (93) 

An = -(7?An - 3(72'An + 3Tr (y^Yu) An + Tr (yJjYn) An 

-2(7?MiYn - e^^^MaYN + 6Tr (y|jAu) Yn + 2Tr (yJjAn) Yn 
+4YnYJjAn + SAnYJjYn + 2YnY;Ae + AnYJ^Ye, (94) 

Ae = -3(7?Ae - 3(72Ae + 3Tr (yJ^Yd) Ae + Tr (y^Ye) Ae 

-6(??MiYe - QglM2YE + 6Tr (yJ,Ad) Ye + 2Tr (y^,Ae) Ye 
+4YeYJ,Ae + 5AeYJ,Ye + 2YeYJ,An + AeYJ,Yn, (95) 

= - Au - SglA^ - ^glA^ + 3Tr (yJ^Yu) Au + Tr (yJ,Yn) Au 

-'^glM^Y^ - GglM^Y^ - ^glM^Y^ + 6Tr (y^Au) Yu + 2Tr (y^An) Yy 
+4YuY{jAu + 5AuY{jYu + 2YuY|,Ad + AuYJ^Yd, (96) 



Ad = --^7?AD-3(72'AD-y(73'AD + 3Tr (yJ^Yd) Ac + Tr (yJ,Ye) Ad 

-^glMiY^ - QglM^Yj, - ^glM^Y,, + 6Tr (yJ,Ad) Yd + 2Tr (y|,Ae) Yd 
+4YDYJ3AD + 5AdYJ,Yd + 2YdY{jAu + AdY{jYu. (97) 
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